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13.9  GHz  RADIOMETER  DESIGN 
V.R.  Reddy 

ABSTRACT 

This  report  concerns  the  design  of  a  13.9  GHz  total  power  radiometer.  A 
large  amount  of  literature  is  available  on  the  subject,  but  this  report 
addresses  a  detailed  design  concept  without  going  into  the  mathematical 
analysis.  All  through  the  discussion,  contrast  Is  brought  forth  between  a 
radar  and  radiometer,  since  the  former  is  a  more  familiar  system  to  most  of 
the  beginners  in  microwave  system  design.  No  stress  is  laid  on  the  theo¬ 
retical  aspects  of  the  radiometer. 

The  system  described  has  an  overall  noise  figure  of  7.67  dB  and  uses  an 
integration  time  of  about  1  sec.  It  is  capable  of  a  resolution  of  2  K  and  the 
slope  of  the  transfer  function  is  set  as  1  mV/1  K  for  convenience..^ 
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1.0  INTRODUCTION 


1.1  Need  for  A  Radiometer 

Radiation  emitted  as  well  as  reflected  by  any  object  is  known  to  be  a 
function  of  various  physical  properties  of  the  object  itself;  hence,  measured 
radiation  data  can  be  used  to  understand  the  properties  of  the  scene  under 
observation.  The  well-established  science  of  measuring  such  radiation  Is 
called  "radiometry"  and  the  device  used  to  carry  out  the  measurements  is  known 
as  a  "radiometer".  This  report  describes  the  systematic  desiqn  and  calibra¬ 
tion  procedure  of  a  13.9  GHz  practical  radiometer  used  to  probe  sea  Ice. 

For  convenience,  the  report  Is  organized  into  six  sections  and  various 
sub-sections.  In  a  broad  sense.  Section  1  describes  the  overall  system. 

Section  2  illustrates  the  practical  design  with  some  analysis  based  on  the 
references.  Section  3  is  used  to  bring  out  some  of  the  methods  of  implement¬ 
ing  the  most-often  required  system  measurements.  Cal ibration  ,of  the  system  is 
treated  in  adequate  detail  In  Section  4,  followed  by  the  description  of 
operation  in  Section  5.  Since  many  of  the  components  used  in  this  system  fell 
short  of  the  requirements,  Section  6  is  used  for  suggestions  to  improve  the 
performance  of  the  system.  Specifications  and  the  measured  data  are  presented 
in  Appendix  A.  Appendix  8  gives  the  FORTRAN  program  listing  used  for 
computing  the  calibration  curves. 

1 .2  Why  Is  It  Different  From  a  Radar  Receiver 

Unlike  a  radar,  a  radiometer  is  a  passive  sensor.  Another  major 
difference  between  a  radar  and  radiometer  lies  in  the  fact  that  the  unwanted 
noise  present  at  the  input  of  the  receiver  antenna  of  a  radar  is  the  signal  of 
interest  in  the  case  of  a  radiometer.  On  paper,  a  radiometer  looks  like  any 
conventional  superheterodyne  receiver,  but  in  practice  the  design  approach  is 


2 


quite  different,  as  will  be  evident  in  a  later  section  on  design.  Basically, 
radiometers  are  divided  into  one  of  two  types:  (a)  total  power  mode  or  (b) 
switched  mode,  the  latter  being  powerful.  The  type  discussed  here  operates  In 
the  total  power  mode,  but  can  be  modified  to  work  in  the  switched  mode. 

1.3  Description  of  a  13.9  GHz  Total  Power  Radiometer 

Figure  1  shows  the  block  diagram  of  the  complete  system.  The  antenna 
used  is  a  standard  gain  horn  connected  to  one  of  the  four  ports  of  a  waveguide 
microwave  switch.  The  receiver  and  noise  diode  are  connected  to  the  other  two 
ports  and  the  last  port  is  left  open.  During  measurements,  the  connections  at 


the  switch  will  be  just  as  shown  in  Figure  1,  i.e.,  the  antenna  is  coupled  to 
the  receiver  and  the  diode  is  coupled  to  the  open  port.  In  this  condition, 
the  switch  is  said  to  be  in  position  1.  During  calibration  (position  2), 
however,  the  antenna  is  coupled  to  the  open  port  and  the  diode  to  the  receiver 
to  deliver  a  known  amount  of  noise.  The  switching  between  positions  1  and  2 
is  done  by  applying  +28  VDC  to  the  switch.  Since  the  isolation  between  the 
two  guides  of  the  switch  is  typically  60  dB,  the  signals  (actually  noise)  in 
the  two  guides  do  not  interact  very  much. 


The  signal  is  amplified  in  the  tunnel  diode  amplifier  and  band-limited  in 
a  narrow-band  filter  with  gooo  (>  50  dB)  stop-band  rejection.  In  practice, 
the  filtering  is  done  in  two  parts,  the  first  one  being  a  lowpass  filter  and 
the  other  a  bandpass  filter.  The  bandpass  filter  is  a  four-cavity  tunable 
filter  and  It  may  also  respond  to  the  harmonics  of  the  signal  (i.e.,  27.8  GHz, 
41.7  GHz,  etc.).  The  lowpass  filter  is  used  to  suppress  these  harmonics.  The 
double-balanced  mixer  used  has  a  waveguide  RF  port  and  10  and  IF  ports  in  SMA 
connectors.  The  local  oscillator  used  for  superheterodyne  mixing  is  a  13.93 
GHz  crystal-stabilized  source.  The  r.iechanism  behind  the  generation  of  the 
13.93  GHz  signal  involves  a  fundamental  oscillator  at  1.393  GHz  which  will  be 
multiplied  by  xlO  to  get  the  required  frequency.  Phase  locking  is  achieved  by 
a  reference  crystal  oscillator  operating  at  100  MHz  which,  in  turn,  will  be 
fed  to  a  harmonic  generator  whose  output  will  be  compared  with  the  1.393  GHz 
fundamental  in  a  phase  detector.  The  error  signal  output  of  the  phase  detec¬ 
tor  is  used  to  tune  the  fundamental  oscillator.  The  isolator  provides  a 
better  match  between  the  oscillator  and  the  mixer.  The  intermediate  frequency 
(IF)  amplifier  is  a  high-gain  discrete  version  operating  at  30  MHz  frequency 
with  a  3  dB  bandwidth  of  only  1.2  MHz.  As  will  become  evident  in  later 
discussion  this  small  bandwidth  is  a  serious  limitation  to  the  feasible  system 
resolution.  Square-law  detection  is  carried  out  by  a  crystal  detector  with 
reasonably  good  sensitivity.  At  this  point,  it  is  interesting  to  note  that 
the  expected  change  in  signal  level  at  the  detector  output  over  the  entire 
dynamic  range  of  temperature  is  just  a  fraction  of  a  millivolt.  This  calls 
for  a  large  video  amplification  to  scale  up  the  signal  to  a  reasonable  value, 
say,  a  few  hundred  millivolts.  The  resulting  signal  is  integrated  for  a 
suitable  period  of  time  (about  a  second)  and  read  out  on  a  digital  voltmeter. 
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1 . 4  How  to  Convert  It  Into  >;  Switching  Type 

To  convert  the  existinq  system  into  a  switching-type  radiometer  the  RF 
switch  must  be  switched  on  and  off  continuously  between  the  antenna  and  the 
reference  load  (the  diode  noise  source  in  this  case).  Usually  a  square  wave 
is  used  for  switching,  thus  spending  50  percent  of  the  time  in  the  desired 
signal  measurement  and  the  other  50  percent  in  calibration.  The  switching 
frequency  is  limited  by  the  RF  switch  response.  Also,  the  RF  detected  signal 
must  be  bandpass-filtered  at  both  the  switching  frequency  and  some  of  Its 
harmonics.  As  a  rule  of  thumb,  at  least  two  harmonics  must  be  included,  and 
the  signal  should  be  synchronously  detected  before  Integration. 

All  of  the  required  microwave  components  were  available  in  the  lab,  but 
not  their  specifications.  Therefore,  every  component  was  evaluated  and  tuned 
wherever  necessary  to  meet  the  required  specifications  as  far  as  possible. 

The  results  of  the  measurements  are  shown  in  Appendix  A. 


2.0  SYSTEM  DESIGN 
2.1  Necessary  Equations 

In  essence,  a  radiometer  must  exhibit  a  linear  transfer  characteristic 
between  the  Input  temperature  and  the  output  voltage.  Since  the  input  power 
is  proportional  to  the  temperature,  it  also  has  a  linear  relationship  between 
the  input  power  and  the  output  voltage.  Mathematically,  the  input  power  to 
the  system,  as  referred  to  the  antenna  terminals  PSy(.,  is  given  by 


sys 


=  k  Tsy$  B 


(1) 


where : 
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k  =  Boltzman's  constant 

B  =  effective  bandwidth  of  the  IF  amplifier  (Hz) 

TSyS  =<  system  temperature  referred  to  the  antenna  terminals  (K) 

In  turn,  the  system  temperature  can  be  expressed  as: 

Tsys  .  nTa  *  (l-n)Tp  *  (L-l)Tp  *  LTrec  (2) 

where : 

n  =  antenna  radiation  efficiency 

L  =  transmission  loss  between  the  antenna  and  the  receiver 
Ta  a  antenna  temperature  of  the  scene  under  observation 
Tp  =  physical  temperature  of  the  system 
Tpec  =  receiver  noise  equivalent  temperature. 

A  detailed  derivation  of  these  equations  can  be  found  in  many  books 
(e.g.,  Ulaby,  Moore  and  Fung,  1981). 

2.2  Set  the  Specifications 

Once  the  input  and  output  relations  are  known,  the  first  thing  towards 
the  system  design  is  to  define  the  following  system  specifications. 

1.  The  dynamic  range  of  the  system,  which  usually  means  the  difference  in 
power  between  the  maximum  signal  (the  fact  that  the  receiver  may  saturate 
at  higher  levels  leading  to  an  undesirable  nonlinear  transfer  function 
sets  the  limits)  and  the  minimum  detectable  signal  (usually  set  by  the 
thermal  noise  floor)  expressed  in  dB.  However,  for  a  radiometer's  system 
temperature  (linearly  related  to  the  power  for  a  fixed  bandwidth  as  long 
as  the  system  is  maintained  at  a  constant  temperature),  it  may  be  more 
realistic  to  express  the  system  dynamic  range. 
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2.  The  required  output  voltage  range  corresponding  to  the  expected  input 
temperature  range,  which  is  the  slope  of  the  system  transfer  function. 

3.  Overall  system  noise  figure,  which  is  the  one  that  sets  the  thermal  noise 
floor  mentioned  in  1  above. 

4.  Bandwidth  of  the  Intermediate  frequency  amplifier. 

5.  Fractional  gain  variation  of  the  overall  system. 

6.  The  dynamic  range  limits  of  the  square-law  detector. 

A  detailed  consideration  of  these  parameters  is  applied  to  the  system 
under  discussion.  Since  the  radiometer  is  a  passive  sensor  and  the  noise 
signals  available  are  very  small,  the  upper  limit  of  the  dynamic  range  Is  not 
a  problem.  In  this  particular  case,  s^nce  the  system  is  expected  to  operate 
at  very  close  ranges  (as  low  as  a  few  meters),  the  lower  limit  also  does  not 
pose  any  real  problem.  In  medical  applications  it  is  more  of  a  problem 

because  the  antenna  will  be  always  In  the  near  field.  In  the  case  of 

radiometers  used  in  astronomy,  where  high  sensitivity  is  of  paramount  impor¬ 
tance,  the  lower  limit  plays  a  critical  role.  For  this  radiometer  a  20K  - 
280K  input  range  is  set  and  the  output  voltages  corresponding  to  this  are 
required  to  fall  between  0  mV  and  260  mV,  which  means  a  slope  of  1  mV/1  K. 

Even  though  the  range  of  temperatures  expected  in  the  case  of  sea  ice  Is  much 

smaller,  the  above  specification  is  set  in  view  of  the  calibration  (which  will 
be  discussed  elsewhere). 

The  measured  overall  system  noi^e  figure  with  the  available  components 
F  =  7.67  dB  =  5.85.  From  this  the  receiver  noise  equivalent  temperature  T 
is  calculated  as: 

Trec  =  (Fn  -  1)  70  «  (5.85  -  1)  273  =  1324. 05K 
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The  gain  variation  aG/G  of  the  complete  front  end  is  assumed  to  be  equal 
to  10  .  In  practice, it  is  somewhat  difficult  to  measure  this  variation  so  a 

setup  with  a  stability  of  at  least  an  order  better  Is  needed.  For  this  dis¬ 
cussion  the  value  assumed  above  is  believed  to  be  reasonable. 

3  dB  bandwidth  of  the  IF  amplifier  available  =  1.2  MHz. 

Dynamic  range  of  the  square-law  detector  *  10  dBm-(-50  dBm)  =  60  dB. 
Other  pertinent  data  for  individual  components  is  shown  in  Appendix  A. 


2.3  Practical  Design 
2.3.1  Front  End 

Using  Equation  (2)  in  the  middle  of  the  dynamic  range,  the  system 
temperature  Tsys  can  be  calculated  as: 

Ta  -  150  K 

Tsys  =  0.7  x  150  +  (1-0,7)273  +  (1.02-1)273  +  1.02  x  1324.05  =  1542.9  K 
Si  mi  1  arly, 

minimum  T$yS  =  1451.9  K 
and 

maximum  Tsys  =  1633.9  K. 


2.3.2  A  Look  at  the  Resolution 

The  minimum  detectable  change  In  antenna  temperature  known  as  radiometric 
resolution  AT,  is  given  by 


AT  = 


y  |_I  +  (— )^  I 
'sys  [Bt  1  rJ  l 


1/2 


(3) 
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The  derivation  of  this  equation  can  be  found  in  Ulaby,  Moore  and  Fung 
[1981]  and  Evans  and  McLelsh  [1977].  The  worst  case  of  resolution  occurs  for 
maximum  Ts^s  and,  hence.  It  is  preferable  to  use  this  value  In  the  above 
formula.  For  different  values  of  integration  time  the  calculated  radiometric 
resolution  is  shown  In  Table  I  as  well  as  in  Figures  2  and  3. 


FIGURE  2:  Radiometric  Resolution  vs  Integration  Time 

for  AG/G  *  .001 


Equation  (3)  shows  that,  for  a  given  system  temperature  and  fixed 
bandwidth,  the  radiometric  resolution  Improves  as  the  square  root  of  the 
Integration  time,  provided  the  fractional  gain  variation  Is  low.  Table  I 
shows  that  the  resolution  saturates  for  AG/G  >  1/Bt.  Hence,  it  Is  very 
Important  that  the  fractional  gain  variation  be  kept  as  low  as  possible  to 
obtain  fine  resolution  for  a  reasonable  value  of  the  time-bandwidth  product. 
Fine-resolution  systems  usually  achieve  that  kind  of  resolution  by  using  large 
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bandwidths  rather  than  large  integration  time.  Also,  since  the  effective 
noise  bandwidth  is  greater  than  the  3  dB  bandwidth  of  the  IF  amplifier,  the 
actual  resolution  in  each  of  the  entries  in  Table  I  will  be  better  (less)  than 
the  values  shown.  For  the  current  system,  since  the  integration  time  is  about 
1  sec,  the  radiometric  resolution  can  be  expected  to  be  approximately  2  K. 

As  can  be  seen  from  equation  (3),  the  radiometric  resolution  AT  Is 
proportional  to  the  system  temperature  when  all  other  quantities  in  the 
equation  are  fixed.  The  Tsys,  In  turn,  depends  upon  the  antenna  temperature, 
which  is  related  to  the  temperature  of  the  scene  under  observation  and  the 
receiver  noise-equivalent  temperature  T„_.  For  fixed  T  the  resolution  AT 

rec  a 

varies  in  proportion  to  the  receiver  noise-equivalent  temperature.  This  Is 
where  the  importance  of  the  low-noise-figure  requirement  for  receivers  (and 
hence  the  use  of  an  expensive  low-noise  amplifier)  comes  into  the  picture. 

In  this  case  for  a  aG/G  of  1CT3  the  effect  of  receiver  noise  figure  on 
resolution  with  integration  time  as  the  parameter  is  calculated.  The  results 
are  shown  in  Table  II  and  also  plotted  in  Figures  4-7. 


TABLE  II 


Receiver  NF 

(dB) 

Radi  ometrl  c 

Resolution  AT  (K) 

t  *  1  msec 

t  =  10  msec 

t  =  100  msec 

t  *  1  sec 

2 

12.90 

4.10 

1.36 

0.60 

3 

16.19 

5.15 

1.71 

0.76 

4 

20.35 

6.47 

2.15 

0.95 

5 

25.58 

8.13 

2.71 

1.20 

6 

32.16 

10.23 

3.40 

•  1.51 

7 

40.45 

12.86 

4.28 

1.90 

8 

50.89 

16.18 

5.38 

2.39 

9 

64.03 

20.36 

6.77 

3.00 

10 

80.57 

25.62 

8.52 

3.78 
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RADIOMETRIC  RESOLUTION  TOO 


RECEIVER  NOISE  FIGURE CdO> 


FIGURE  4:  Radiometric  Resolution  vs  Receiver  Noise  Flqure 

for  t  -  1  mSEC 


FIGURE  5:  Radiometric  Resolution  vs  Receiver  Noise  Flqure  for 

t  -  10  mSEC 


Table  III  and  Figure  8  show  the  dependence  of  radiometric  resolution  AT 
on  the  noise  bandwidth,  provided  all  other  quantities  are  held  constant.  The 
plot  is  calculated  from  equation  (3)  using  maximum  T$yS,  t  of  1  sec,  and  aG/G 
-  10*2  an(j  io~3. 


TABLE  III 


Bandwidth  B 

at  (k) 

AT  (K) 

(MHz) 

with  AG/G  =  10-3 

with  AG/G  =  lO-2 

- 1 - - 

T5755 

2 

1.91 

15.61 

3 

1.80 

15.60 

4 

1.74 

5 

1.71 

6 

1.68 

7 

1.66 

8 

1.64 

9 

1.63 

10 

1.61 

This  Indicates  that  no  matter  how  much  bandwidth  and  "t"  you  use,  a  gain 
variation  of  10~2  is  totally  unacceptable.  In  fact,  for  a  1  K  resolution  and 
a  feasible  Integration  time  of  1  sec,  AG/G  needs  to  be  much  better  than  10”3 
(maybe  10-4).  The  other  alternative  is  to  use  a  switching  mode  radiometer. 

2.3.3  IF  Section  and  Detection 

The  next  step  In  the  design  process  is  to  calculate  the  required  gain  for 
the  complete  front-end  system  (i.e.,  up  to  the  square-law  detector).  The  cri¬ 
terion,  as  mentioned  earlier,  is  to  see  that  the  signal  (noise)  power  corres¬ 
ponding  to  the  mid-point  of  the  temperature  dynamic  range  coincides  with  the 


% 
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FIGURE  8:  Radiometric  Resolution  vs  Bandwidth 
for  AG/6  «  .001 

mid-point  of  the  square-law  detector  dynamic  range.  It  should  be  noted  that, 
If  no  other  component  prior  to  the  detector  is  driven  In  to  saturation,  the 
dynamic  range  Is  going  to  be  limited  by  the  detector.  However,  in  radiometers 
this  situation  seems  almost  Impossible  as  long  as  the  above-mentioned  cri¬ 
terion  Is  satisfied.  In  the  present  case,  the  ratio  of  the  maximum  to  the 
minimum  system  temperature  Is  nuch  less  than  2,  indicating  that  we  need  less 
than  3  dB  of  the  detector's  available  dynamic  range.  Nevertheless,  it  Is  good 
practice  to  position  the  signals  at  the  center  of  the  response  of  the  square- 
law  detector.  In  this  case  an  antenna  temperature  of  130  K,  corresponding  to 
a  system  temperature  of  843.88  K,  must  produce  a  -20  dBm  (10*^  W)  signal  at 
the  input  of  the  square-law  detector.  The  receiver  gain  required  to 
satisfy  this  condition  would  be  calculated  as  follows: 
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V 


(4) 


where : 


rf  "  K  T"  B  Fn 
sys  n 


6r^  =  receiver  power  gain  (without  antenna)  in  dB 

=  power  output  of  the  IF  amplifier  corresponding  to  the  T$yS. 


rf  1.3805xl0‘23  x  1542.9  x  1.2xl0+6x  5.85 


=  6.69x10  '  =  78.25  dB 


Note  that  If  the  actual  effective  noise  bandwidth,  rather  than  the  3  dB 
IF  bandwidth,  is  used  in  the  expression,  Grf  may  decrease  about  a  dB  or  so. 
For  this  system,  the  tunnel  diode  amplifier  has  a  gain  of  16  dB  and  the 
single-side  band  conversion  loss  of  the  mixer  is  21.2  dB.  The  pre-selection 
filter  Insertion  loss  is  0.3  dB  and  the  waveguide  loss  between  the  amplifier 
and  the  switch  is  0.1  dB.  Therefore,  the  required  IF  amplifier  gain  Gif  1s 
given  by: 


G1f  ■  G  f  -  16  +  21.2  +  0.3  +  0.1  =  83.85  dB. 


A  crystal  detector  of  sensitivity,  r,  500  mV/mW  is  used  for  square-law 
detection.  The  sensitivity  of  such  a  detector  is  expressed  in  millivolts/ 
millwatts  because  the  Input  power  is  linearly  related  to  the  output  voltage, 
not  to  the  output  power.  From  this  information  the  minimum  and  maximum  values 
of  the  voltages  available  at  the  input  of  the  video  amplifier  are  calculated 


vin  •  r  x  K  *  Tsys  x  8  x  Fn  *  Grf 


Minimum  V.  =  500x1. 3805xl0"23xl451. 9x1. 2xl0+6x5. 85x6. 69x10*''  »  4.707  mV 
i  n 
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Maximum  V.  =  500x1. 3808xl0_23xl633. 9x1. 2xlO+6x5. 85x6. 69xl0+7  =  5.297  mV 
i  n 


2.3.4  Back  End 

Notice  that  the  change  in  the  detector  output  voltage  is  just  0.32  mV  for 
a  change  of  260  K  at  the  Input  of  the  antenna.  This  calls  for  a  large  video 
amplification  to  scale  up  the  small  voltage  change  to  a  sufficiently  large 
change.  Thus  the  video  amplifier  gain  Gv  is  computed  as: 


desired  output  slope 
available  input  slope 


1  mV/K 

'5.297  mV  -  4.707  mV)/ (280  K 


=  440.678 


Gv  =  20  log  (440.68)  *  52.88  dB. 


Thus  at  the  output  of  the  video  amplifier  the  voltage  levels  are: 


20  K  +  2.074  V 

280  K  ♦  2.334  V 


Finally,  by  providing  an  offset  voltage  of  2.074  V,  we  get 


0.0  mV 


280  K  -*■  260  mV 
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Another  requirement  is  that  the  video  amplifier  must  provide  proper 
terminating  resistance  for  the  detector.  In  this  case  it  should  be  approx¬ 
imately  1  K.  Also,  because  of  the  dc  component  present  (the  actual  signal  is 
not  dc,  but  a  very  low  frequency  which  will  be  further  integrated)  at  the 
output  of  the  detector,  the  final  offset  might  be  somewhat  different  from  the 
calculated  1.089  V  value.  Thus  It  Is  advisable  to  provide  the  offset  after 
integration  rather  than  before. 

As  mentioned  above,  all  of  the  front  end  components  were  available  in  the 
lab.  The  measured  data  for  each  of  these  components  are  listed  in  Appendix 
A.  The  complete  back  end  of  the  radiometer  is  built  as  a  single  module.  The 
circuit  diagram  is  shown  in  Figure  9.  It  consists  of  a  high  qain  (approxi¬ 
mately  60  dB)  video  amplifier  IC-1,  lowpass  filter  IC-2,  and  an  offset  circuit 
IC-3.  An  instrumentation  amplifier,  rather  than  discrete  operational  ampli¬ 
fiers,  is  chosen  for  video  amplification.  It  not  only  reduces  the  chip  count 
but  also  eliminates  any  possibility  of  small  oscillations  associated  with  high 
gain  amplification  (with  high  gain  even  a  small  positive  feedback  is  enough  to 
initiate  oscillations).  These  oscillations  may  sometimes  be  very  close  to  the 
desired  signal  frequency,  making  it  impossible  to  eliminate  even  by  lowpass 
filtering.  Since  an  instrumentation  amplifier  is  Internally  compensated,  all 
of  these  problems  will  automatically  be  taken  care  of.  An  additional  advan¬ 
tage  of  the  Instrumentation  amplifier  Is  the  fact  that  the  IC  has  an  on-chip 
common  point  (operational  amplifiers  do  not  have  this)  making  It  less  susceo- 
tible  to  spurious  pickups.  As  shown  in  Figure  9,  potentiometer  R1  is  used  to 
adjust  the  offset  of  the  instrumentation  amplifier  and  R2  to  trim  its  gain. 

The  input  is  terminated  in  1  K  ohms  to  provide  matching  to  the  output  of  the 
square  law  detector.  The  integrator  or  lowpass  filter  is  built  with  a  BIFET 
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FIGURE  9:  Video  Amplifier  Circuit  Diagram 


operational  amplifier  (LF  355)  and  has  a  time  constant  of  RC  -  1.9+®  x  0.47“^ 

»  0.893  seconds.  It  is  a  two-pole  Butterworth  filter.  Finally,  the  offset 
circuit  is  realized  with  a  BIFET  as  a  differential  unity  gain  amplifier. 
Potentiometer  R3  is  used  to  adjust  the  final  offset  voltage.  Thus,  the  output 
of  this  circuit  is  the  required  integrated  signal  and  read  on  a  digital 
voltmeter. 

All  the  modules  are  mounted  on  a  2-mm-thick  aluminum  plate.  The  layout 
of  the  various  modules  and  interconnections  is  shown  schematically  in  Figure 
10.  Also  shown  are  the  dc  power  supplies  required  for  different  modules. 


3.0  SYSTEM  FCASUREHENTS 

Two  of  the  most  Important  specifications  of  any  microwave  reciever  system 
are  power  gain  and  noise  figure.  Each  component  of  the  system  will  affect,  to 
some  extent,  one  or  both  of  these  quantities.  Hence,  proper  measurement  of 
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FIGURE  10:  Layout  of  the  Various  Modules 


these  quantities  for  each  of  the  system  components  (to  make  required  adjust¬ 
ments  to  the  system)  and,  finally,  the  system  (to  evaluate  the  system  perform¬ 
ance)  is  warranted.  Usually,  at  microwave  frequencies  three  types  of  measure¬ 
ment  techniques  are  available:  (1)  manual;  (2)  semi-automatic;  (3)  automatic. 

For  single-frequency  measurements,  manual  techniques  are  quite  adequate; 
semi-automatic  techniques  are  better  suited  for  swept-frequency  measure¬ 
ments.  If  high  accuracy  is  required,  automatic  techniques  are  used.  The  last 
two  are,  naturally,  much  faster  and  often  much  easier  to  use.  In  this  case. 

It  Is  a  single-frequency  measurement  and  manual  techniques  are  employed. 

3.1  Power  Gain  or  Loss 

In  one  of  the  simplest  setups  (Figure  11)  a  siqnal  source  and  a  power 
meter  are  used.  As  usual,  some  Interconnecting  cables  and  adaptors  are 
required.  At  the  frequency  of  interest  (13.9  GHz)  the  si^.ial  input  to  the 
device  under  test  (OUT)  In  dBm  Is  measured  by  bypassing  the  OUT  as  shown. 
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Next,  the  output  of  the  OUT  Is  noted  down  in  dBm.  The  difference  in  the  two 
readings  gives  the  gain  or  loss  in  dB.  There  are  many  drawbacks  to  this 
setup,  but  the  one  of  most  concern  is  the  mismatch  between  any  two  components 
of  the  setup  --  source  and  OUT,  OUT  and  power  meter.  The  reflections  caused 
by  these  mismatches  result  In  Inaccurate  readings.  Between  the  source  and  the 
DUT  an  isolator  could  be  used  (VHF  through  microwaves)  but  between  the  OUT  and 
the  power  meter  it  may  not  be  possible  because,  If  the  OUT  is  the  complete 
receiver  ar.d  the  output  is  measured  at  the  IF  stage,  no  isolators  may  be 
available.  Quite  often,  however,  some  amount  of  attenuation  is  mandatory 
between  the  DUT  and  the  power  meter  in  order  to  protect  the  power  meter  from 
being  damaged  by  the  larqe  outputs  of  the  OUT  (Figure  12).  This  automatically 
takes  care  of  the  mismatch  problem  to  some  extent.  Instead  of  an  isolator,  an 
attenuator  could  be  used  as  long  as  the  signal  level  Is  within  the  specifi¬ 
cations  of  the  DUT.  In  addition,  a  directional  coupler  of  known  coupling  is 
used  to  measure  the  Input  signal  to  the  OUT  without  disconnecting  it  from  the 
setup.  The  amount  of  the  input  attenuation  required  is  decided  by  the  level 
of  matching  desired.  For  example,  if  there  is  a  VSWR  of  2.0:1  between  the 
source  and  the  DUT  and  a  VSWR  of  1.2:1  Is  desirable,  the  needed  attenuation  is 
half  the  dB  difference  betwen  the  two  return  losses  —  one  corresponding  to 
each  VSWR.  So,  magnitude  of  reflection  coefficient 

without  the  attenuator  =  (2-l)/(2+l)  *  0.33 

return  loss  =  -20  log  (0.33)  =  9.54  dB 

magnitude  of  reflection  coefficient 

with  the  attenuator  =  (1.2-l)/( 1.2+1)  =  0.09 

return  loss  =  -20  log  (0.09)  =  20.9  dB. 

Therefore  the  minimum  required  attenuation  =  (20.9-9. 54 ) / 2  =  5.68  dB. 
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FIGURE  11:  Basic  Setup  for  Power  Gain  Measurement 
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FIGURE  12:  Improved  Setup  for  Power  Gain  Measurement 


A  rigorous  analysis  (Including  the  multiple  reflections)  yields  a 
slightly  larger  value.  Sometimes  a  variable  attenuator  Is  required  to  adjust 
the  Input  power  level  to  the  OUT  in  which  case  the  mismatch  problem  Is 
automatically  taken  care  of. 
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3.2  Noise  Figure 

The  method  described  here  Is  known  as  the  Y-factor  method  and  this  Is,  by 
far,  the  most  often  used.  One  way  of  Implementation  is  shown  in  Figure  13. 
Noise  generated  by  a  diode  noise  source  of  known  excess-noise  ratio  (ENR )  Is 
fed  to  the  OUT.  An  Isolator  of  known  insertion  loss  Is  used  to  reduce  the 
mismatch.  Once  again,  an  attenuator  with  low  VSWR  could  be  used  Instead.  The 
IF  amplifier  (30  MHz  or  60  MHz)  Is  not  required  If  the  DUT  Is  the  complete 
receiver  and  DUT  output  is  the  IF  amplifier  output  of  the  receiver.  The 
variable  attenuator  is  used  to  adjust  the  power  level  and  also  to  reduce  the 
mismatch.  The  power  meter  may  have  to  be  very  sensitive  (the  one  used  here 
has  a  range  of  -60  dBm  to  -20  dBm). 

When  the  diode  Is  off  (cold  condition  -  0  VUC)  the  output  power  Is  noted 
In  dBm  as  N1  and  the  diode  Is  switched  on  (hot  condition  -  +28  VDC).  The 
output  Is  noted  again  In  dBm  as  N2.  Now  the  noise  figure  In  dB  is  calculated 
as  [Laverghetta,  1981]: 

NF  =  ENR  (dB)  -  Insertion  loss  of  the  Isolator  - 
10  log  [Antilog  (N2/N1/10)  -  l]  dB 


FIGURE  13:  Setup  Noise  Figure  Measurement 
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If  an  attenuator  is  used  in  place  of  the  isolator  its  attenuation  must  be 
used  In  place  of  the  Isolator  Insertion  loss  in  the  above  equation.  The 
contribution  of  the  external  IF  amplifier  could  be  removed  from  the  OUT  noise 
figure  if  the  IF  amplifier  noise  figure  and  the  OUT  gain  are  predetermined. 

4.0  CALIBRATION 
4.1  Different  Methods 

Like  any  other  linear  system,  radiometer  calibration  Involves  determining 
the  relation  between  the  output  and  input,  in  other  words  the  system  transfer 
function.  In  the  case  of  a  radiometer,  the  output  parameter  Is  the  voltage 
measured  after  Integration  and  the  input  parameter  Is  the  system  tempera¬ 
ture.  Since  the  system  Is  linear,  at  least  two  coordinates  on  the  transfer 

characteristic  must  be  found  In  order  to  complete  the  straight  line.  The 

« 

source  of  the  Input  must  be  stable  and  its  antenna  temperature  be  known  quite 
accurately.  This  is  very  Important  because  the  whole  system's  performance 
depends  on  this.  It  may  also  be  a  good  idea  to  choose  the  two  coordinates  so 
that  'rh  one  is  at  a  different  edge  of  the  transfer  characteristic.  There 
are  many  alternatives  that  can  satisfy  these  requirements.  Some  of  them  are: 

1.  Clear  (scatter  free)  sky 

2.  Eccosorb 

3.  Still  water 


4.  Noise  diode. 


4.2  How  to  0o  in  Practice 


The  reference  sources  1,  2  and  4  are  quite  dependable  calibration 
sources.  For  a  given  frequency,  the  antenna  temperature  of  the  sky  is  a  func¬ 
tion  of  the  angle  of  observation,  atmospheric  quantities  (pressure,  tempera¬ 
ture,  water  vapor  density  profiles),  and  the  presence  of  scattering  media  like 
clouds.  For  a  clear  sky  condition  all  of  the  parameters  above  can  either  be 
measured  or  estimated  with  good  accuracy.  Therefore,  this  is  one  of  the  most 
commonly  used  reference  loads  In  radiometry.  Depending  on  the  various  para¬ 
meters  discussed  above,  the  sky  temperature  may  be  anywhere  between  2  K  and  30 
K.  A  FORTRAN  program  to  calculate  this  temperature  with  different  Inputs  is 
available  on  the  Honeywell.  For  reference,  the  listing  is  given  in  Appendix 
B.  The  graph  In  Figure  14  is  obtained  for  standard  atmospheric  profiles  using 
this  computer  program.  To  obtain  the  first  coordinate,  point  the  radiometer 
antenna  towards  the  sky  at  some  convenient  angle  (small  angles  from  zenith  are 
preferable  because  the  sky  temperature  doesn't  vary  too  much).  Make  sure 
there  are  no  clouds  In  the  sky  and  note  the  output  voltage.  The  other  coor¬ 
dinate  may  be  obtained  by  blocking  the  antenna  aperture  with  eccosorb,  a 
microwave  absorbing  material  which  Is  close  to  perfect  ( 1 . e . ,  100%)  at  these 
frequencies;  It  Is  well  known  that  a  perfect  absorber  exhibits  an  emlsslvity 
of  unity.  This  means  the  antenna  temperature  will  be  the  same  as  Its  physical 
temperature  as  long  as  the  absorber  Is  kept  close  to  the  antenna.  The  phy¬ 
sical  temperature  of  the  absorber  can  be  measured  accurately  by  embedding  a 
thermocouple  (a  metallic  heterojunction  which  develops  a  voltage  across  the 
junction  proportional  to  the  junction  temperature)  in  the  absorber.  The 
output  voltage  measured  in  this  condition  gives  the  second  coordinate.  The 
straight  line  joining  these  two  coordinates  (Figure  15)  can  be  used  to  relate 
any  output  voltage  to  the  corresponding  system  temperature. 


Still  water  Is  another  useful  reference  source  if  the  properties  of  the 
water  (dielectric  constant,  physical  temperature,  foam  coverage,  etc.)  are 
known  accurately.  Also  required  Is  the  sky  temperature  (discussed  above). 
The  computer  program  mentioned  before  can  also  be  used  to  calculate  the  an¬ 
tenna  temperature  of  still  water  [Ulaby,  Moore  and  Fung,  1984],  Note  this 
temperature  Is  also  polarization  dependent.  Figure  16  shows  the  plots 
obtained  using  the  computer  program  for  vertical  as  well  as  horizontal  polar 
Izatlons.  Once  again,  small  angles  measured  from  nadir  seem  to  be  better 
suited  for  calibration. 


NAOIR  ANCUE (DEGREES) 


FIGURE  16:  Hater  Temperature  vs  Nadir  Angle  at  13.9  GHz 


Reference  4  has  some  advantages  over  1  and  3  in  the  sense  that  it  is  not 
atmospheric  dependent.  In  this  case  a  diode  capable  of  producing  very  wide 
band  noise  is  turned  "on"  and  "off"  with  a  +28  VDC  supply.  If  the  noise  power 
generated  by  this  diode  in  one  condition  (on  or  off)  is  known,  the  noise  power 
in  the  other  condition  can  be  calculated  provided  the  excess  noise  ratio  of 
the  diode  is  known.  In  turn,  it  can  be  related  to  temperature  and,  hence,  can 
be  used  for  calibration.  This  type  of  reference  is  Ideal  for  switching  types 
of  radiometers.  Also,  it  may  be  noted  that  by  calibrating  the  instrument  at 
regular  intervals,  the  effect  of  the  system  gain  variations  on  the  desired 
output  can  be  minimized. 


5.0  OPERATION 

Operation  of  this  system  is  fairly  simple.  Like  many  other  systems,  this 
radiometer  needs  some  warm-up.  Therefore,  switch  on  the  system  and  keep  on 
standby  for  at  least  10  minutes.  This  assures  a  stable  local  oscillator  out¬ 
put  within  the  specified  limits  and  also  minimizes  drifts  in  other  circuits. 

If  the  RF  switch  is  not  already  In  position  2,  apply  +28  VDC  to  the  appro¬ 
priate  terminals  (see  Appendix  A)  to  change  from  position  1  to  position  2. 
Point  the  antenna  towards  the  sky  at  some  angle  and  observe  the  integrated 
video  voltage  output.  If  It  is  drifting  too  much,  leave  some  more  time  for 
warm-up  and  see  whether  there  Is  any  improvement  In  the  output  stabiity.  If 
it  is  not  within  the  expected  acceptable  limits,  reduce  the  IF  amplifier  gain 
slightly  by  turning  the  gain  control  potentiometer  in  a  counter-clockwise 
direction.  A  full  revolution  of  the  potentiometer  changes  the  gain  by  as  much 
as  20  dB,  therefore,  a  fraction  of  the  full  revolution  should  be  attempted  at 
a  time.  Note  that  the  output  voltage  will  change  and  so  does  the  offset 
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requirement.  Now  place  the  eccosorb  in  front  of  the  antenna  and  note  the 
output  voltage.  From  this  data,  calculate  the  output  slope,  then  trim  the 
video  amplifier  gain  (increase)  slightly  to  get  the  desired  slope.  Repeat 
this  until  the  expected  stability  is  reached.  At  this  point,  adjust  the 
offset  using  R3.  The  receiver  is  now  ready  for  measurements. 

For  satisfactory  operation,  calibration  must  be  performed  as  often  as 
possible  at  regular  Intervals. 

6.0  COMMENTS  ON  THE  SYSTEM 

6.1  Antenna 

The  horn  antenna  used  in  this  system  exhibits  very  high  sidelobe  levels, 
leading  to  an  antenna  temperature  Ta  consisting  of  the  contribution  from  the 
mainlobe  as  well  as  from  the  sidelobes.  In  simple  terms,  this  means  the 
received  signal  is  not  only  from  the  observation  area  (covered  by  the  main 
beam),  but  also  from  the  larqe  surrounding  area  (covered  by  many  sidelobes  to 
different  extents),  which  is  clearly  an  undesirable  effect.  In  the  case  of  an 
Imaging  system,  this  situation  becomes  even  more  undesirable  because  the  side¬ 
lobes  have  the  effect  of  overlapping  the  scene  under  observation.  Kw,;-  1,  it 
i*  reasonable  to  use  low  sidelobe  antennas  like  a  reflector  antenna,  an  array 
antenna  or  (better)  a  corrugated  horn. 

6.2  Low  Noise  Amplifier  (LNA) 

Recall  that  the  low  noise  amplifier  sets  t.he  lower  limit  of  the  overall 
system  dynamic  range  because  this  amplifier  is  designed  to  have  a  lower  noise 
figure  (i.e.,  lower  thermal  noise  floor)  than  the  rest  of  the  system  that 
follows.  Also,  it  is  well  known  that  the  NF  contribution  of  the  subsystem 


29 


(excluding  LNA)  to  the  complete  system  (including  LNA)  will  be  less  by  a 
factor  of  the  amplifier  gain  (actually,  a  little  better  than  this  -  for  exact 
expressions  refer  to  Ulaby,  Moore  and  Fung  [1981]).  It  is  for  this  reason  the 
low  noise  amplifier,  usually,  is  the  first  component  in  any  receiver.  An 
exception  to  this  Is  that  in  high  power  radars  (transmitting  power,  not  the 
rec  '  r4  echo  signal)  using  a  single  antenna  for  transmission  and  reception,  a 
pro  _  device  called  a  "limiter"  (sometimes  two  of  them:  one  a  tube  type 
and  the  other  a  solid-state  type)  will  be  used  before  the  amplifier  to  protect 

4 

the  sensitive  receiver  from  large  leakage  power  coupled  from  the  transmit¬ 
ter.  In  radiometers,  there  Is  no  need  for  any  such  protection,  but  a  switch 
is  usually  required  before  the  amplifier  for  calibration  purposes.  Also,  the 
use  of  such  LNA  is  always  advantageous  in  radar  if  properly  selected,  whereas 
in  radiometers  it  has  demerits  as  well . 

Merits: 

1.  A  smaller  noise  figure,  lower  thermal  noise  floor,  thus  a  better 
sensitivity. 

2.  Allows  the  use  of  slightly  inferior  components  In  the  rest  of  the  front 
end. 

Demerits: 

1.  If  not  properly  selected,  the  upper  limit  of  the  system  dynamic  range  may 
be  degraded. 

2.  The  foremost  disadvantage  is  the  large  gain  variations  associated  with 
these. 

3.  They  are  more  expensive  than  any  other  component  in  the  system. 

4.  Above  millimeter  frequency  range,  no  solid-state  LNAs  are  currently 
available  in  the  market. 
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In  light  of  the  above  discussion,  for  this  particular  system  it  may  be 
reasonable  to  eliminate  the  LNA,  if  a  good  quality  mixer  (low  conversion  loss 
and  high  LO  to  RF  Isolation,  among  other  things)  is  available  for  use. 

6.3  Hi  xer-Preampl 1 f 1 er 

The  mixer  currently  used  has  very  bad  high  single  sideband  (SSB) 
conversion  loss  and  very  low  LO  to  RF  port  Isolation.  The  undesirable  effect 
of  the  low  isolation  Is  not  only  that  some  of  the  signal  is  lost  in  multiple 
reflections,  but  also  that  the  variations  in  the  local  oscillator  noise 
contribution  will  appear  as  a  signal  at  the  receiver  output  leading  to 
erroneous  data.  Hence,  the  system  performance  can  be  improved  by  using  a 
better  quality  mixer. 

Currently  no  preamplifier  Is  used  in  this  system.  But  a  mixer  integrated 
with  a  preamplifier  capable  of  a  combined  gain  (preamp  gain-mixer  conversion 
loss)  of  about  30  to  40  dB  is  desirable. 


6.4  IF  Amplifier 


Assuming  that  all  of  the  other  parameters  remain  the  same,  the  radio- 
metric  resolution  could  be  Improved  (reduced)  by  increasing  the  bandwidth  (see 
Section  2.0).  For  a  fixed-time-bandwidth  product  the  Integration  time  can  be 
reduced  by  increasing  the  bandwidth  thus  achieving  not  only  the  same  resol¬ 
ution,  but  also  higher  output  data  rates.  Because  of  this,  radiometers 
(unlike  radars)  always  operate  with  a  large  3  dB  bandwidth.  In  this  case  a  3 


dB  bandwidth  of  8  MHz  is  not  only  reasonable  but  also  quite  practical,  as  many 
manufacturers  supply  30  MHz  IF  amplifiers  with  as  much  as  10  MHz  bandwidths. 


6.5  Square-Law  Detector 


Proper  selection  of  this  device  is  critical  for  radiometers.  Some  of  the 
important  features  to  be  considered  are:  dynamic  range,  linearity  of  the 
transfer  characteri Stic ,  sensitivity,  and  stability.  The  first  three  can  be 
satisfied  by  most  of  the  detectors  on  the  market,  whereas  the  last  one  limits 
the  choice  of  detectors.  Crystal  detectors  have  good  sensitivity  but.  In 
terms  of  stability,  they  are  not  the  best.  Tunnel  diode  detectors  are  known 
to  demonstrate  almost  an  order  better  stability  compared  to  crystal  detectors, 
but  with  reduced  sensitivity.  While  the  reduction  in  sensitivity  can  be 
compensated  for  elsewhere  In  the  design,  the  stability  cannot;  therefore.  It 
is  much  better  to  use  a  tunnel  diode  detector  for  square-law  detection. 

6.6  Video  Amplifier 

The  fact  that  the  temperature  coefficient  for  various  passive  components 

■ 

(resistors  and  capacitors)  used  in  this  system  can  cause  remarkable  changes  in 
the  resolution  of  the  system  is  important.  Hence,  it  is  reasonable  to  use 
very  low  tolerance  (may  be  1%)  components. 


7.0  SUMMARY 

A  detailed  and  systematic  design  procedure  for  a  practical  13.9  GHz 
radiometer  has  been  presented.  Various  design  tradeoffs  were  discussed  and 
the  effect  of  various  parameters  on  the  feasible  resolution  of  the  total  power 
radiometer  was  shown.  Required  modifications  were  spelled  out  for  converting 
this  unit  into  a  switching-type  radiometer.  To  make  the  task  of  trouble¬ 
shooting  easier,  some  of  the  important  measurements  were  also  discussed.  A 
simple  calibration  procedure  was  used  to  illustrate  the  principles.  Many  of 
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the  difficulties  with  the  components 
clear  and,  in  turn,  some  suggestions 
receiver  performance. 


used  in  the  current  system  were  made 
were  prescribed  to  improve  the  overall 
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APPENDIX  A 


Standard  Sain  Horn 
Model  number 
Nominal  gain 
Nominal  bandwidth 

Directivity 
Radiation  efficiency 


12-12 
24.7  dB 
9°  E-plane 
10°  H-plane 
26.23  dB 

70.3%  at  13.9  GHz 


Figure  A-l  shows  the  dimensions  of  the  antenna. 


FI6URE  A.l:  Antenna  Dimensions 


Waveguide  Switch 

Model  number 

82152-33000200 

Insertion  loss 

0.1  dB 

Isolation 

60  dB 

VSWR 

1.1:1 

Actuator  voltage 

+28  V 

Current 

1A  max. 

Wiring: 

Inside  the  switch 

A  *  Black 

Outside  the  switch 

B  ■  Orange 

C  *  White 

A  »  Black 

B  *  Red 

C  »  White 

Figure  A-2  shows  the  switch  positions  and  the  connections  required  to  change 
from  one  position  to  the  other. 


FIGURE  A.2:  Waveguide  Switch  Details 


Tunnel  Diode  Amplifier 


Model  number 
3  dB  bandwidth 
Nominal  gain 
Noise  figure 
VSWR 

Power  output  at  1  dB  GCP 
Supply  voltage 
Current 


T8503 

13.4  -  14.4  GHz 
16  dB 
4.85  dB 
1.1:1 
-4.0  dBm 
-15  VOC 
15  mA 


No  adjustments  are  required. 


Bandpass  Waveguide  Filter 

Bandwidth 

Insertion  loss 

Stop  band  attenuation 


50  MHz 
0.3  dB 

>  50  dB  (100  MHz  away  from  3  dB  point) 


Tuned  for  maximum  receiver  gain 


Crystal  Oscillator 

Model  number  SSX-0225  (El) 

Frequency  13.929  GHz 

Power  output  6.2  dBm  (measured  at  the  output  of  the 

isolator) 


Supply  voltage 


-20.0  VOC 


Model  number 

1210-2 

DSB  noise  figure 

17. 

2  dB 

LO-RF  Isolation 

8.2 

dB 

VSWR  at  the  RF  port 

1.6 

:1  at  13 

.9 

GHz 

IF  Amplifier 

Model  number 

EVT30LL37 

Center  frequency 

30 

MHz 

3  *.rf  bandwidth 

1.2 

MHz 

Nilse  figure 

2.0 

dB 

ndxlmum  power  gain 

82 

dB 

Saturated  power  output 

+22 

dBm 

Supply  voltage 

+15 

VDC  at 

170 

1  mA 

-15 

VDC  at 

30 

mZ 

Crystal  Detector 
Model  number 
Sensitivity 


Dynamic  range 


HP  423B 
500  mV/mW 
60  dB 


Video  Amplifier 


Nominal  gain 


Integrator  time  constant 


Input  impedance 


60  dB 
.893  sec 


IK  ohm 


APPENDIX  B 
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ft 


ai?40 

01220 
31260 
01270 
31230 
0 1 2-70  20 

3I330C 
21210 

21220 
•  i  •*«•** 

*j  *  j 

2 1  243  <77 

21250  23U 
3122.3  231 
01270 

oizcec 

31270 
01400 
31410 
01420 
01420b 
01440 
01420 
314*00 
01470 
01400 
01490 
01300C 
01310 
01320 
01320 
0134OC 
01330 
01360 
01370 
01380 
01390b 
01600 
01610 
01620 
01620b 
01640 
01630b 
01660 
01670b 
31630 
01690b 
01700  10 

01710C 
B1720 
01720 

01740 

31730 

O1760C 

01770 

017S0C 

01790 

01800 

01010 

01320 

01320C 

01840 

01830C 


TAU  -  KA  *  DELTA 

TAU02  -  TAU0Z  *  3. 3* (TAU  *  TAUPl 
TAUP  -  TAU 

IF(Z.LT.H*.0001.OR.  2.ST.H-.  3001  >  TAUOH-TAU0Z 
CONTINUE 

SUN  -  SUM  ♦  KA*T.EXP<-1.  »rAlj0Z»SEC)  «C 
COMPUTE  Tdowo  AND  raky 
TDN  -  SUM  •  DELTA.  »  SEC 
TGKY  *  TDN  *  2.7  •  EXP ( - 1 . »TAUOZ »2E2) 

tuaTm7Auoh*t0* ( i .  -TAUJi-n  *  <  i  .  csr>  *ts-  < 1 .  -auoi>  »gsp»tsky 
PE  TURN 

F3PMAT  flX,'  ALT  1.02  :■  S  :~C  k  WEICHT') 

FORMAT . IX, F4. 1 . 4F3.3) 

END 

SUBROUTINE  H2OVAP0R 
SUBROUTINE  H20VAPQR ( F , P , T , R , KH2Q ) 

REAL  K)  120 

G  *  2. 33*F/1013. *(200. /T) **0.626*1 1. +0. 018*R*T/P> 

KH20  *  2.  ♦F»F»R*(20O. /T>  »*l.3*G*  (200.  /T»£XP(  -64*.  /T)  /  t  (494.  4 
-F»F)*-*2*4. *F*F*G*G!  *  1.2C-0O) 

RETURN 

ENO 

SUBROUTINE  OXYGEN 
SUBROUTINE  OXYGEN  (F,P,T.R.K02> 

REAL  K02,  FN (23, 2) ,  YN(20,2> 

COMMON  /OLOCKl/FN,  YN 

INITIALIZE  THE  CONSTANTS 
GN  -  1.10  *  P/1012,  •  (200. /T) **0.83 
GB  -  0.49  *  P/1013.  *  (200.  /T)  *-*0. 89 
SUM  »  0. 

PERFORM  SUMMATION  OVER  N 

I  -  a 

DO  10  N-l ,29,2 
I  -  I  ♦  l 

PH IN  *  4.6C-03  ♦  200. /T  •  (2.*N  ♦  1.)  * 

EXP(-6.B9G-03  *  N  *  (N*l. ) *200./T7 

DNP  -  SORT (N  •  (2. *N  *  3. )/((N  *  l.)*(2.*N  ♦  !.))> 

DNN  »  SORT ( (N  *  l.)*(2.*N  -  t.)/(N*(2.*N  *  l.))> 

GN1  -  <GN*ONP*DNP  *  P*(F-FN(1 , I) ) *YN(I,1) 7/ 

< (F-FN ( 1,1*) **2  *  GN*GN) 

GN2  6  (GN* DNP* DNP  *  P* (-1 .  *F-FN (1,1)) *YN (1,1))/ 
<(-l.*F-FN<I,J))**2  *  GN*GN) 

GN3  *  (GN* DNN* DNN  *  P* (F-FN ( I , 2) ) «YN (I , 2) ) / 

( (F-FN (1,2)) **2  *  SN*GN) 

GN4  »  (GN* DNN* DNN  *  P* (-1 .  *F-FN ( 1 ,2) >  *YN < 1 ,2) ) / 

( (-1 . *F-FN( 1,2) ) **2  *  3N*GN) 

SUM  -  SUM  *  PH IN  *  (GN1  ♦  GN2  ♦  GN3  *  3N4) 

COMPUTE  F‘  b  K02 

FPRIME  -  0,7  •  3B/(F*F  *  GB*GB)  ♦  SUM 

K02  -  1.61E-02  •  F*F  •  FPRIME  »P/1013.  *  <3@0./T>*»2 

RETURN 

EKr- 

yilBROUTINE  ATMOS 

SUBROUTING  ATMOS ( Z , P, T ,R, PO . TO . RO) 

COMPUTE  THE  TEMPERATURE  AT  ALTITUDE  Z 
Til  -  T0  -  6.3*11. 

T  -  TU 

IF  (Z.LT.ll.)  T  -  T0  -  6.3*1 
IF  (Z. 8T.20. )  T  -  Til  ♦  Z  -  20. 

COMPUTE  THE  PRESSURE  AT  ALTITUDE  Z 
P  -  P«  /  EXP1Z/7. 7) 

COMPUTE  THE  WATER-VAPOR  DENSITY  AT  ALTITUDE  Z 


Best  Available  Copy 


1 

aao ibc 

MAIN  CALLING  ROUTING 

00020 

REAL  FN(20,2> ,  YN(20,2) 

00030 

INTEGER  ■'OL 

00040 

LOMPLEX  EW 

aooso 

COMMON  /DLQCX1/FN,  YN 

OOOoO 

REAL  KH20 .  K.02 

•30070 

DATA  {FN.1,1)  .1-1, 20) /36. 2640,53. 4466, 39. 3710, 60. 4340, 

300001. 

61. 1306.61.2002,62.  4 112. 62.  9790, 63.  3633, 64.  1270,64.6727, 

oocro*. 

63.2241 ,63.  '647,60.  . *323,66.2367, 67.2^74,67.  7007,63.4300, 

aoiooi 

63.7601 ,67. 4CC 7/  • 

aaiio 

DATA  (FN  v  1 ,2)  ,1*1,201/1. 1C.  7303,62.  4063,60.  3061  ,37.  1642, 

oo i  rot. 

33.3237,37.  j  123,36.  7632,36. 3634, S3.  7033,33.  2214 ,34. 671 1 , 

oa : 30*. 

54.  1. 30, 33 .  j757, 5. .3 .6 C ,  52 .  —  422 ,32.3212,31.33-.O,^0.  > C .  . , 

30143.* 

S3.  473... 47.  76 13/ 

00  ISC 

DATA  iYN( I , l> , I»l ,20)  /  4.  3 1. 4.  74, 3.  52, 1.06,. 333,-1. <33, 

00  l  63* 

*2.23, -3.32,  -4. 32,  -3. 26,  o.  13,  -6.77 , -7. 74, -8. 61 , -9. 1 1 , 

00:70!. 

-10.3,  7.37,  -  13.2,  ••7.87,-23.0/ 

00100 

DATA  (YN.t ,2)  . 1*1 ,20) /-O. 214, -3.70, -3. 92,-2. 68, - 1 . 13, 

00: 70S. 

.344,1.63,2.04,3. 71 ,4.93,3.04,6.76,7.53,3.47,7.01 ,10.3, 

00  ZOO!. 

7.26, 13.3,7.31 ,26. 4/ 

00210c 

SCALE  THE  YN  '  S  PROPERLY 

00223 

DO  3  1*1,20 

oozra 

YN<I,l>  *  7H ::,1>  •  1 E -34 

00240 

3 

YN( 1,2)  »  YN»:,2)  *  IE-04 

0025OC 

INPUT  SEQUENCE 

00260 

triUE  -  02 

00270 

WR I TE  ( 06 , 20 1 ) 

00200 

READ (03. 200)  IX 

30270 

IF  < IX-3)  ,99,99 

CC2O0 

IF  (IX-l)  6,6, 

30310 

WRITE (06, 202) 

00320 

GO  TO  99 

30330 

6 

WR I TE ( 06 . 203 ) 

00340 

READ (03,200)  P0,T0,RB 

00330 

IF  (PO)  ,  ,  7 

00360 

WRITE (06, 204)  "PRESSURE* 

00370 

GO  TO  6 

00300 

7 

IF  (T0)  ,,a 

00390 

WRITE (06 , 204)  “TEMPERATURE* 

00400 

GO  TO  6 

00410 

3 

IF  (RO)  ,  ,  9 

00420 

WRITE (06, 204)  "WATER  VAPOR* 

00430 

GC  TO  6 

00440 

9 

WR I TE ( 06 , 203  > 

00450 

00460 

r?;r:;r  F  Best  Available  Copy 

00470 

WR  ITE  (06 , 204)  "FREQUENCY" 

00480 

GO  TO  9 

00490 

11 

IF  (F-100.)  10,10, 

00300 

WRITE (06, 206) 

00310 

GO  TO  9 

00320 

10 

WRITE (86 ,210) 

00330 

READ (33, 200)  IX 

00340 

IF  (IX-l)  12,12, 

00330 

PRINT, -ENTER  POLARI ZATIONi l-H, 2-V ’ 

00360 

REAO  20, POL 

00370 

20 

FORMAT (V) 

00580 

PRINT, 'ENTER  REAL  V  I  MAG  PARTS  OF  WATER  DIEL  CONST- 

00390 

READ  20, EW 

00600 

PRINT, ’ENTER  RADIOMETER  ALTITUDE  <K«)  A DOVE  WATER  LEVEL' 

00610 

REAO  20, H 

00620 

12 

WRITE (06, 207) 

00630 

READ (03, 200)  A NO 

00640 

IF  (ANG)  ,13,13 

80630 

WRITE (06,204)  "ANGLE" 

ata&70  i3 
306G0 
30670 
00700  14 

30710 
00730 
30733 
00740  79 

00733  130 
3073 1  % 
007331. 
'307603. 
337704 
007714 
307S0  330 
00790  301 
000004 
303104 
OO030  333 
00330  303 
008404 
300304 
00860  304 
30870  303 
00880  306 
008904 
30900  307 
00710  308 
309304 
0O9T0  210 
009404 
00930 
0O960C 
00970 
00980 
00990 
01  TOO 
01010C 
01030 
01030 
01040 
01030 
01060 
01070 
01080 
01090 
01100 
01  110 
01130C 
01120 
01140 
01  130 
01160 
01170 
01180 
01  190 
ai300C 
01310 
01320 
01320' 


12 

IF  (ANC-QO.)  14,14. 

WK:TEC0&,30O) 

CO  TO  13 

CALL  TEMP  (P.rOL,34,H,ANG,P0, 70 ,R0, TSKY , TWAT, TAU02 ,  ITTLC) 
WR1TE(QL.  1(30)  F  , PCI. , 34 ,11, FO  ,  TO , .90 ,  T0X7  ,  TWAT . TALOZ 
WPrTE(  IFIL2,  130)  P  ,  POL  ,  £W ,  H  , P3  ,  T0  ,R0  ,  TSKY ,  TWAT ,  T ALIO 2 
CC  TO  .2 
crop 

format ;  is.  rr.ca  -  .rs.  i ,  ghs  • , 

/  ,  IX,  Pol  jri:  IH^l  ,'/*;)  -*',12,'  Uat^r  dial  con*t-»  '  ,  3F6. 3 , 

/  ,  l::  ,  '  A1  tx  Cude*  ‘  ,  r'3.  2,  Km  , 

/,l  A, 'Pa  *  '  ,  F6 .  1 ,  Jitur  To  »',F3.1,'  X  Ro  •  '  ,  F3.  3 ,  ’  g/.ft2  '  , 

/,1 X,  T ®ky  -  ',F9.4,'  K  T»*tar  -  ,F7.4,'  X', 

/,  IX, ‘OPTICAL  DEPTH  (Tju)  -',Fa.3,'  Np  '  > 

FORMAT (V) 

FORMAT (IX,/, IX, 'ENTER  1  FOR  STANDARD  ATMOSPHERE ‘ l X , 

3  FOR  USER  DEFINED  ATMOSPHERE ' ,/,  1 X , 

2  TO  EXIT  PROGRAM') 

FORMAT ( IX, 'MODI FY  PROGRAM  FOR  USE  WITH  NEW  ATMOSPHERE  MODEL') 
FORMAT (IX,  ENTER  SURFACE  PRESSURE  (*6*rt) ,  SURFACE  TEMPERATURE', 
'  <0 agraa*  Xal vtn) , '/ , IX, '  AND  SURFACE  WATER  VAPOR  DENS IT/  ', 

'  (g/imG)  ' ) 

FORMATdX,/,  'INVALID',  All,  ENTRY',/,  IX,  TRY  AGAIN') 

FORMAT ilX,' ENTER  FREQUENCY  <GHx)'> 

FORMATdX, 'FREQUENCY  SHOULD  NOT  EXCEED  100  GHj’,/,lX, 

•—  —  PLEASE  TRY  AGAIN') 

FORMATdX, 'ENTER  ANGLE  FROM  ZENITH  (Dagra**)') 

FORMATdX,  ANGLE  IS  NOT  *3  EXCEED  80  Dagraa* ' ,  / ,  IX , 

' -  PLEASE  TRY  AGAIN') 

FORMATdX,  "ENTER  l  FOR  Taky  COMPUTATION"  ,/,  IX  , 

*  2  FOR  Twatar  CALCULATION  AT  A  GIVEN  ALTITUDE") 

END 

SUBROUTINE  TEMP 

SUBROUTINE  TEMP  ( F ,  POL ,  EW  ,  W ,  ANG ,  P0  ,  T0  ,  R0  ,  TSJCY  ,  TWAT , TAU0Z , IFILE) 
REAL  X02 ,  X!  '20 ,  KA 
INTEGER  POL 
COMPLEX  £W 

INITIALIZE  THE  CONST  ANTS 
PI  -  4. »ATAN ' • . ) 

sec  -  i./con  /*hg*p:/i80. ) 

N  -  2200 
Z0*>  22. 

DELTA  -  Z0/FLOAT(N> 

SUM  -  0. 

TAUP  ■  0. 

TAU0Z  *  0. 

WR I TE  d  F I  L£ ,  30®  > 

CALL  3AMA(POL,EW,H,ANG,3SP) 

BEGIN  INTEGRATION 
DO  20  1-1  ,N-H 
C  *  l. 

IF  (I.EO. l.OR. I.EO.N*!)  C  -  3.3 

Z  -  (I-l. ) »OELTA 

CALL  ATMOS(Z,P,T,R,P0,T0,R0) 

CALL  H20VAP0R ( F , P , T , R , XH20 ) 

CALL  OXYGENIF ,P , T,R ,K02> 

CONVERT  FROM  dB/k»  TO  Nap ar/ka 
K02  -  X02  •  0.23026 
XH20  -  FH20  *  0.22026 
XA  •  X02  ♦  XH2D 


Best  Available  Copy 


01960 
0197a 
01980 
01870 
01700 
01710 
01711 
01720 
01720 
01740 
Q1730C 
01760 
01770 
017S0C 
01990  10 
02000 
00010  20 
02020 


ft  »  90  /  EXP (2/2. ) 

return 

END 

CULP  OUTING  GAMA(POL,EN,il,  ANG.GSP) 

COMPLEX  EW , A , 9 

INTCCER  POL 

AN-  AN8»2. 14137/ 103. 

A— C2QRT  (CW-2IN  (Af.  >  >3  IN  (AN) 

3-CMPLX (COO  (AN> ,q. ) 

IF IrGL. EQ. 2)  oCTO  10 

horizontal  polarization 
GSP««2AC0 ( (C- A) :  (2*Ai  > > **Z 
CCTO  20 

VERTICAL  POLARIZATION 
CONTINUE 

GSF- (CAES  (  (  (EW*D)  -  A)  /  (  <£W»E)  *A)  >  )  »-»2 
RETURN 
END 


Best  Available  Copy 


